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We study the couphngs between collective vibrations such as the isovector giant dipole and 



CN ■ isoscalar giant quadrupole resonances in tin isotopes in the framework of the time-dependent 

> ; 

. Hartree-Fock theory with a Skyrme energy density functional. These couplings are a source of 

I anharmonicity in the multiphonon spectrum. In particular, the residual interaction is known to 

I couple the isovector giant dipole resonance with the isoscalar giant quadrupole resonance built on 

■ top of it, inducing a nonlinear evolution of the quadrupole moment after a dipole boost. This 
coupling also affects the dipole motion in a nucleus with a static or dynamical deformation induced 

^ I by a quadrupole constraint or boost respectively. Three methods associated with these different 



manifestations of the coupling are proposed to extract the corresponding matrix elements of the 
residual interaction. Numerical applications of the different methods to ^^^Sn are in good agree- 
ment with each other. Finally, several tin isotopes are considered to investigate the role of isospin 
and mass number on this coupling. A simple 1/A dependence of the residual matrix elements 
is found with no noticeable contribution from the isospin. This result is interpreted within the 
Goldhaber- Teller model. 
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I. INTRODUCTION 



A particular interest in strongly interacting systems is their ability to present disorder or 
chaos, and, in the same excitation energy range, well-organized motion. Atomic nuclei are 
known to show both behaviors [1]. In particular, they exhibit a large variety of collective 
vibrations, also called giant resonances (GRs), with excitation energy usually above the 
particle emission threshold [2]. The GRs are associated with anomalously large cross sections 
in some nuclear reactions. 

Baldwin and Klaiber observed the isovector giant dipole resonance (GDR) in photofission 
of uranium nuclei [3], interpreted as a vibration of neutrons against protons [4]. This GDR 
was investigated with several probes [2] and was also observed on top of highly excited 
states, e.g., in hot nuclei [5]. The survival of ordered motion in hot nuclei, i.e., in a chaotic 
environment, is one of the most striking phenomena in nuclear physics. Other kinds of GR 
have been discovered, such as the isoscalar giant quadrupole resonance (GQR) associated 
with an oscillation of the shape between a prolate and an oblate deformation [6], and the 
isoscalar giant monopole resonance (GMR) corresponding to a breathing mode [7-9] . 

The GRs are usually associated with the first phonon of a small-amplitude harmonic 
motion. However, the proof of their vibrational nature came with the observation of their 
two- and three-phonon states [10-12]. Multiphonon studies also provided a good test of the 
harmonic picture. In particular, anharmonicity was found in an abnormally large excitation 
probability of these states, indicating that different phonon states couple because of the 
residual interaction [13, 14]. Microscopic investigations, such as the random phase approxi- 
mation (RPA) together with boson mapping techniques [15] and the nonlinear response to an 
external field in the time-dependent Hartree-Fock (TDHF) theory [16, 17] showed, indeed, 
that strong couplings between GMR, GQR and GDR occur. In particular, a GMR or a GQR 
(resp. a GMR) can be excited on top of a GDR (resp. a GQR), leading to couplings between 
one- and two-phonon states. As a consequence, GRs cannot be described in a purely har- 
monic picture. Anharmonicities were also found to affect pygmy dipole resonances, though 
depending on the choice of the nuclear functional [18]. 

The goal of the present work is to get a deeper insight into the couplings between various 
GRs, which represents a first step toward understanding complexity and disorder in nuclei 
at high excitation energies. As an example, we focus on the coupling between isovector 
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dipole and isoscalar quadrupole vibrations. A clear link between the linear dipole motion 
on a deformed state and the quadratic response of the quadrupole moment to an external 
dipole excitation (investigated in Ref. [16]) is made. The TDHF theory is used to compute 
the residual interaction coupling the one-phonon state of the GDR to the two-phonon state 
with a GQR built on top of the GDR. Applications to spherical tin isotopes are performed 
to investigate the role of the isospin degree of freedom and of the total number of nucleons 
on the coupling. 

We present a schematic model describing couplings between GRs and their effect on one- 
body observables in Sec. II. The TDHF formahsm and its apphcation to nuclear vibrations 
are discussed in Sec. 111. Numerical details on the 3-dimensional TDHF code are also given. 
A detailed investigation of the couplings in ^^^Sn is presented in Sec. IV, together with a 
more systematic analysis in tin isotopes. Finally, we conclude in Sec. V. 



II. A SCHEMATIC MODEL FOR GR COUPLING 



Let us illustrate the effect of couplings between vibrational modes within a simple 
schematic model introduced in Ref. [16]. We consider the Hamiltonian 

H^Ho + V (1) 

where corresponds to the harmonic (HF+RPA) part and the residual interaction V 
couples collective modes. Eigenstates of Hq are one- and two-phonon states lu) and 
with eigenenergies E^, — Eg + hu!,^ and E^^ — Eo + hu!i, + hu!ij_, respectively, where cUn^i, denote 
the collective frequencies and Eg is the ground-state energy. In the following, h is omitted 
in the notation. Only the coupling between the two states |z/) and |z//i) is considered here. 
The associated matrix element of the residual interaction is noted = (z/|y|z/yu). Such 
couplings between one- and two-phonon states have been proven to be the most important 
one in nuclei [15]. Using first-order perturbation theory, the eigenvalues of H are those of 
Hq with eigenstates 

\V) - e^\uii) (2) 

and 

[IT/Z) fti (3) 

where Su. — —- 
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The couplings are expected to affect the evolutions of expectation values of one-body 
observables such as the multipole moments Qi,{t) = {Qu){t). We investigate below three 
different manifestations of the couplings on these evolutions. They will be used in the next 
section to compute from TDHF calculations in the case of coupling between giant dipole 
and quadrupole resonances. 

A. Quadratic response 

The effect of couplings in the quadratic response has been introduced in Ref. [16]. High- 
hghts on the main steps are given here. At initial time, the ground state |0) of the system 
is excited by a boost with the one-body operator Q^, 

|*(0)) = exp(-iA;,Q.)|0). (4) 

Developing the exponential up to second order in the boost intensity ki, and considering an 
evolution under the Hamiltonian defined in Eq. (1), the state at time t reads at first order 
in 

ex.p{-iEot) 

where Qi, — (z/|(5,y|0) is the transition amplitude that we assume to be real. 

The expectation value of the one-body observable used in the boost exhibits oscillations. 
Indeed, in case of no static deformation in the ground state, we have 

Q^{t) = -2k^qlsm{uja) + 0{kl). (6) 

In particular, its amplitude increases linearly with the boost intensity in the small amplitude 
regime. In addition to this linear response, the coupling induces an oscillation of the other 
collective mode Q^: 

Q^{t) ^ 2klqlq,^ [cos^i) - 1] (7) 

where we have assumed = {fJ'\Q^\0) = {fiiylQuliy) ■ This oscillation is then quadratic in 
and provides a first method to compute the residual interaction v^, assuming the fact that 
a nonhnear theory, such as TDHF, is used to follow the expectation values of the one-body 
observables. We finally note that Qiy{t) and have different frequencies and start in 

phase quadrature. 
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|0) - ik^q^e-'^"' (|I7) - e^e-'''^'\i7Jl)) 



(5) 



B. Linear response in an external static field 



It is interesting to note that the couphng may also manifest itself in the linear response 
to the boost (4) if an external static field is added to the Hamiltonian (1) 

^(A)=^(0)+Ag^. (8) 

We choose A small enough to induce a linear static deformation defined as 

Ql{X) = {0{X)\Q,\0{X))^x(?^^^ ^, (9) 

where the ground state |0(A)) of ^(A) contains a contribution of the one-phonon state 

The external potential modifies linearly the eigenenergies of the Hamiltonian and the fre- 
quency of the linear response to a boost (4) on |0(A)) follows 

providing another direct way to extract the matrix element of the residual interaction. We 
emphasize the fact that, here, the nonlinear response is not invoked and a RPA code allowing 
static deformation in the ground state would be sufficient to compute such couplings. 



C. Response to two simultaneous excitations 

We showed two manifestations of the coupling (i) in the quadratic response and (ii) in 
the linear response under a static constraint. Let us now introduce a third one where the 
response Qu{t) is studied after a simultaneous double boost of and Q^: 

1^(0)) = e-'''''^>'e-'^''^''\0). (12) 

The Qn term modifies the response of Eq. (6) with an additional term 

= 4:k^ki^qlqi^-^[l - cos{uJi^t)]cos{ujj). (13) 
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It is convenient to write this evolution with the form 

— sin uj^t — P cos uj^t + ^cos{uJi^ + uj^)t + ^cos{ujfi — uj^)t (14) 

where Q^, = —2kyql and {3 = 2k^q^Vij^/uj^. In fact, we can show that x{t) is a solution of the 
differential equation 



X 

+ 



1 - 2/3^sina;.t 



,2 



X + 13— COS ujfjt = Q (15) 



if one keeps only the first-order terms in The two first terms of the left-hand side 

are equivalent to a Mathieu's equation. It is not surprising because the latter has been 

shown to qualitatively reproduce the preequilibrium dipole motion coupled to collective 

shape vibrations of the system in N/Z asymmetric fusions [19, 20]. 

We see in Eq. (14) that the effect of the coupling produces vibrations at frequencies \iOi,±. 

Lijfj\. By analogy to the standard response function related to the strength distribution [21] 

we introduce the coupling response function 

—1 

Rl{uj) = — — / dt cos{u;t)AQ,{t) (16) 

defined for cj > 0. The latter can be used to investigate the coupling because it is linearly 
proportional to v^: 

2 

Rtiuj) « [-2(5(0;, - a;) + S{uj^ + uj^ - uj) + S{\u, - - u)] . (17) 

Equation (17) provides then a third way to extract v^. 

Let us finally note that the contribution to the coupling response function at cUi, and those 
at jo;,^ ± have opposite signs in Eq. (17) and that the integral of the coupling response 
function is zero. It is interesting to note that this property is still valid at all order in k^, and 
kf^. To show it, let us recall that, in our schematic model, Qu\0) = QuW) and Ql\0) = ql\0), 
which implies e"*'^"*^" |0) = c,|0) — iSulu), where c^, = cos{ki,qi,) and = sm{ki,q^). The 
response Qu{t) following the double boost of Eq. (12) then becomes 

V 

Qu{t) = -2cvS„qv s\n.{ijj^,t) + Aci^c^SuS^qu— [1 - cos(a;^t)] cos{u^t). (18) 

We see that Eq. (14) is still valid if one replaces Q^, by Q'^, — —2ci,Si,qi, and (3 hy j3' — 
2Cfj_Sfj,Vfj_/uj,j_. Then, the a;-dependance of Eq. (17) is unchanged. As a consequence, the 
cancellation of the integral of the couphng response function defined in Eq. (16) is not 
limited to the small-amplitude regime. 
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III. THE TIME-DEPENDENT HARTREE-FOCK APPROACH 



A. Applications to nuclear vibrations 

Coherent motion of fermions such as collective vibrations in nuclei can be modeled by 
time-dependent mean-field approaches like the TDHF theory proposed by Dirac [22] . Indeed, 
in its hnearized version, TDHF is equivalent to the RPA which is the basic tool to understand 
the collective vibrations in terms of independent phonons. 

As we saw in the previous section, giant resonance properties can be investigated by 
studying the response of the system to an external (collective) one-body field. In particular, 
time evolution of one-body (collective) observables, which can be computed using mean-field 
approximations, contain the necessary information to investigate the couplings between col- 
lective modes. Indeed, TDHF takes into account the effects of the residual interaction if 
the considered phenomenon can be observed in the time evolution of a one-body observable. 
In particular, the nonlinear response in TDHF contains the couplings between one- and 
two-phonon states coming from the 3-particlc 1-hole and 1-particle 3-hole residual interac- 
tion [16]. In that sense, it goes beyond the RPA, which is a harmonic picture and contains 
only 1-particle 1-hole residual interaction. 

In its unrestricted form (i.e., with no constraint on spatial symmetry), TDHF autho- 
rizes all possible spatial forms of the nucleon wave functions, which is crucial because of 
both shell the effects and the wave dynamics. In addition, Landau spreading and evapora- 
tion damping arc well accounted for [23]. However, it does not incorporate the dissipation 
from two-body mechanisms [24-26]. Inclusion of pairing correlations is possible within the 
time-dependent Hartree-Fock-Bogolyubov theory [27] , but realistic apphcations in three di- 
mensions are not yet achieved. Extension to theories going beyond the one-body limit such 
as extended TDHF [26], second RPA [28, 29], time-dependent density matrix theory [30-32] 
or stochastic one-body transport theory [33] should be considered for realistic description of 
giant resonance properties [34]. 

Application of TDHF to nuclear dynamics has been possible thanks to the Skyrme- 
type effective interaction [35, 36]. Early reahstic TDHF codes have been applied to study 
collective vibrations in nuclei with simplified Skyrme interactions [37]. Recent increase of 
computational power allowed realistic TDHF description of giant resonances in 3 dimensions 
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with full Skyrme energy density functional (EDF) [16, 38-41]. In particular, TDHF has been 
used to investigate nonlinear effects in nuclear vibrations [16, 42]. 



B. Formalism 



The TDHF equation can be written as a Liouville-Von Neumann equation 



(19) 



where p is the one-body density matrix of an independent particles state with elements 

p{rsq, r's'q') = fl Vi{rsq)cp*{r's'q'), (20) 
1=1 

where A is the number of nucleons. The sum runs over all occupied single-particle wave 
functions (pi and r, s, and q denote the nucleon position, spin, and isospin respectively. 
The Hartree-Fock single-particle Hamiltonian h[p] is related to the EDF, denoted by E[p\, 
through 



C. Numerical details 

In this work, the TDHF equation (19) is solved iteratively in time on a spatial grid 
with a plane of symmetry using the tdhfSd code built by P. Bonche and coworkers [43] 
with the SLy4 parametrization of the Skyrme EDF [44]. The latter has been constrained 
on the pure neutron matter equation of state to improve the description of exotic nuclei. 
For instance, together with a density-dependent zero-range pairing interaction, it allows a 
somewhat good reproduction of the isotopic shifts between proton and neutron mean-square 
radii in lead isotopes. It also improves the description of the isotopic evolution of the binding 
energies [44] . Defining the neutron drip line as the isotope for which the chemical potential 
vanishes, it is estimated to be around ^'''^Sn in the tin isotopic chain [45], though it might 
depend on the choice of the pairing interaction and on the inclusion of beyond mean-field 
correlations. 

Good convergences of the quadrupole and dipole moment evolution is ensured with a 
lattice spacing A r = 0.6 fm and a time step A t — 5 x 10~^^ s. The size of the half-box 




(21) 
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where the single particle wave functions are evolved is 80 x 80 x 40 in mesh size unit A r, 
unless otherwise specified. 



IV. RESULTS 

Let us now investigate the coupling between isovector dipole and isoscalar quadrupole 
vibrations in tin isotopes in the framework of the theoretical model presented in Sec. II where 
\v) = \D) and |//) = \Q) denote a GDR and a GQR phonon respectively. The isovector 
dipole moment is defined as 

NZ 

Qj,^^(Zn-Zp) (22) 

where Z„ (resp. Zp) measures the neutron (resp. proton) average position on the z axis. 
The isoscalar quadrupole moment reads 

QQ-^ti^^^-£^-y^)■ (23) 

Their expectation value evolutions are computed using the tdhfSd code after different 
initial conditions as described below. 

A. Nonlinear quadrupole motion induced by a dipole boost 

We first investigate the quadratic response presented in Sec. II A in the ^^^Sn nucleus. 
Figure 1(a) shows the early time evolution of the dipole moment after a dipole boost accord- 
ing to Eq. (4) in the small-amplitude regime. The dipole moment follows a — sin function 
as indicated by Eq. (6). Extracting the frequency from the first minimum of D{t) leads 
to a GDR energy oi lod ~ 15.2 MeV. This value is shghtly lower than the maximum of 
the experimental GDR peak energy E^nax — 16.1(7) MeV [46]. The same analysis in ^^'^Sn, 
in which almost all the dipole strength is located around the GDR energy, gives a value 
of uj\) ?a 15.3 McV which is in good agreement with experimental data where a peak 
energy of E^^'^ = 15.4 MeV has been obtained [47]. Note that, the extraction method of 
cui, from the first extremum of Qu{t) is, in first approximation, comparable to the ratio of 
the second over the first energy weighted moments of the strength function 1712/ mi [16]. 

We see in figure 1(b) that an oscillation of the quadrupole moment is induced by the dipole 
boost. According to the theoretical model presented in Sec. II, this is a manifestation of the 
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FIG. 1: Time evolution of the dipole (a) and quadrupole (b) moments in ^^^Sn after a dipole boost 
with an intensity ko = 0.01 fm~^. 

residual interaction of Eq. (1) coupling the dipole and quadrupole vibrations. In particular, 
Qqit) starts in phase quadrature with Qoit) and oscillates with a smaller frequency. These 
observations are in qualitative agreement with the quadratic response in Eq. (7). 

To get a deeper insight into this coupling, we have computed the TDHF response for 
several dipole boost velocities ko- The first extrema of the dipole and quadrupole moments 
are reported in figure 2(a) and (b) respectively. Whereas the dipole amplitude is indeed linear 
mko StS expected from equation (6), indicating that these calculations are performed in the 
small-amplitude regime, the induced quadrupole motion is quadratic in koi in agreement 
with Eq. (7). 

To obtain a quantitative estimate of the coupling, we first extract the transition amplitude 
from a linear extrapolation of Q™" at /cd — s> in Fig. 2(a). According to Eq. (6), we get 

~ 6.98 fm. The same analysis with a quadrupole boost in the linear regime gives a 
transition amphtude qq ^ 61.4 fm^ and a GQR energy of loq fa 13.0 MeV. Note that the 

I9n (^^^Sn) 

same analysis in "Sn gives a GQR energy of a;^ ' ^ 13.3 MeV, in excellent agreement 
with the experimental value Eq^'j^ — 13.24 ± 0.13 MeV [48]. These quantities, together 
with a quadratic extrapolation of the quadrupole maximum at fc^ ^ in Fig. 2(b), give, 
according to Eq. (7), a matrix element of the residual interaction Vq^ ~ —0.61 MeV. 
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FIG. 2: Circles: First minimum and maximum of the dipole (a) and quadrupole (b) moment 
evolution, respectively, in ^^^Sn as a function of the dipole boost intensity ko- Dashed lines: 
linear and quadratic extrapolations at fco — ^ of the dipole (a) and quadrupole (b) amplitudes, 
respectively. 

B. Dipole motion in a nucleus with a static quadrupole constraint 



The formalism developed in Sec. II B, where the linear response is investigated in an 
external potential, cannot be directly applied to study the coupling between the dipole and 
quadrupole modes. The reason is that the external potential —XQq with the definition of 
Eq. (23) is not bound from below and its use in constrained HF calculations would lead to 
unphysical results. It is then necessary to consider another external potential such as 



where 



Qm — 



E 

i=l 



£.2 



(24) 



(25) 



is the monopole moment and kx = \/b/2 if A > and — \/5 if A < 0. The expression (24) 
then reads 

' ^2 if A > 0, 



3 —A y 

V 167r ^1 _£2 _ ^2 if A < 0. 



(26) 
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FIG. 3: (a) Static quadrupole moment from HF calculation (circles) under a quadrupole+monopole 
constraint (see text) as a function of the Lagrange parameter A in ^^^Sn. (b) TDHF energy of 
the GDR (circles) from the first minimum of the dipole moment after a dipole boost along the 
deformation axis with an intensity = 0.01 fm~^. Dashed lines : linear extrapolations at A — >^ 0* 
of the quadrupole moment (a) and GDR energy (b). 

Such an external field allows one to explore all quadrupole deformations from oblate 
(A > 0) to prolate (A < 0) shapes as shown in figure 3(a) where the ground state quadrupole 
deformation Qq of the constrained HF solution is plotted as a function of the Lagrange 
parameter A. The quadrupole deformation is clearly linear in this perturbative regime and 



its slope at the origin is 



A=0 



-1260.4 fm^ MeV-^ 

As discussed in Sec. II B, such a static deformation is expected to change the dipole 
frequency as compared to that of the GDR excited on the spherical ground state. In fact, the 
frequency of a dipole oscillation along the main quadrupole axis decreases (resp. increases) 
with a prolate (resp. oblate) deformation. This is indeed what we observe in figure 3(b) 
where the energy of the GDR is plotted as a function of A. Note that, according to Eq. (11), 
this is consistent with the negative sign of the ratio vq/qq obtained in Sec. IV A. 

We also observe in Fig. 3(b) that the evolution of this energy is linear both for A > and 
A < 0, but the slopes are different in these two regimes. This is attributed to the presence of 
the monopole moment in the constraint (24). Indeed, the monopole vibration is also coupled 
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to the dipole mode by a matrix element vm of the residual interaction [15, 16]. According 
to Eq. (11), the dipole energy is expected to be modified as 



A compression of the nucleus increases the dipole frequency, which implies vm/qm < 0. 
Because Xkx > for all A, the monopole and quadrupole moments have an opposite effect 
on for A < and act in the same direction for A > 0. This is indeed what we observe 
in Fig. 3(b) where the effect of the constraint almost cancels on the prolate side, whereas it 
strongly increases the GDR energy in the oblate one. 

Finally, starting from Eq. (27), it is possible to isolate the coupling matrix element 
between the dipole and quadrupole modes 



Using the data extracted from Fig. 3 and the value of qq obtained in Sec. IV A, we get 

(2) 

Vq ~ —0.56 McV. This result is in reasonable agreement with the one obtained with the 
quadratic response. 

C. Response to a dipole+quadrupole boost 

A third manifestation of the coupling between dipole and quadrupole motions occurs 
when both a dipole and a quadrupole boost are performed at initial time. We showed in 
Sec. II C that, in such a case, the dipole motion is affected by the quadrupole vibration. Such 
effect is not present in the linear response theory because the modifications arc proportional 
to kokq. As can be seen in Eq. (6), there is no other quadratic term because the next-order 
terms affecting the dipole motion are in k^ and kg. 

The basic tool to study the effect of the coupling on the dipole motion is the couphng 
response function defined in Eq. (16). In principle, its calculation implies to follow the dipole 
moment over an infinite time. However, we use a filtering procedure to avoid numerical 
artifacts coming from the interaction of the nucleus with reflected nucleon wave functions 
because of the hard box boundary conditions [49]. We perform the calculations over 3000 
iterations in time (450 fm/c). The dipole moment is multiplied by a filtering function 




(27) 




(28) 




This procedure induces an additional width of 
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FIG. 4: Coupling response function of the dipole moment after a dipole+quadrupole boost with 
intensities ko = 0.01 fm~^ and kq = 0.001 fm~^ respectively. 

~ 2 MeV. According to Eq. (17), this additional width is sufficiently small for the present 
discussion because the modes in the coupling response function are located at lod—^q ~ 2.2, 
liJd 15.2, and a;^) + a;Q 28.2 MeV. However, the low-energy part of the spectrum, i.e., in 
the region of the loj^ — loq peak, is dependent on the choice of the filtering function within 
these numerical conditions. We checked with other filtering functions, e.g., a cosine instead 
of a Gaussian function, to confirm that the higher part of the spectrum (above 10 MeV) 
is not affected. In addition, the filtering function does not change the fact that the total 
integral of the coupling response function vanishes (see Sec. II C). The latter was found to 
be a solid numerical property of this function. Finally, we checked the convergence of the 
results presented in this section by comparing with calculations performed in a bigger box 
of 120 X 120 X 60 in mesh size unit A r. 

Figure 4 shows the coupling response function for the dipole motion following a 
quadrupole+dipole boost. We checked that, in the small-amplitude limit, the coupling 
response function is indeed independent of kq and ko- As expected from Eq. (17), two 
peaks are present in this energy range at a;^ and lod^ ^jJq with opposite signs. Moreover, 
the integral of the positive peak at a;^) is directly related to the coupling as 



With the coefficients calculated in Sec. IV A, we obtain a couphng v'^^^ fa —0.68 MeV of the 
same order of magnitude than with the two previous methods. 




(29) 
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Let us finally note that in the case of more complicated vibrations, e.g., the oscillations 
exhibit several frequencies, the coupling response function can be used for a more detailed 
investigation of the coupling. Indeed, it allows an analysis of the coupling effect at each 
energy whereas the two previous methods give only access to a weighted sum of the matrix 
elements of the residual interaction associated to each excited mode [16]. 

D. Evolution of the coupling with isospin and mass 

We now repeat the study of the linear quadrupole motion induced by a dipole boost, 
described in Sec. IV A, to the tin isotopic chain. The choice of this method to investigate 
more systematically the coupling between dipole and quadrupole vibrations is motivated 
by its rather low computational time as compared to the two other methods. Our goal is 
to understand the evolution of vq as a function of the isospin. To avoid any ambiguity 
coming from possible static deformation in the ground states, we focus on some of the tin 
isotopes that are spherical at the HF level: ioo.i06,ii4,i2o,i32,i40gj^ These isotopes allow for 
an investigation of the coupling from the proton-rich to the neutron-rich side. 

Let us first investigate the linear response to a quadrupole boost (Eq. (4)) to compute 
the energies uoq and transition amplitudes qq from the first minimum of the quadrupole 
moment (see Eq. (6)). These quantities are plotted in Fig. 5 as a function of the number 
of nucleons. The GQR energy is known to be proportional to A~^/^ [2]. This is compatible 
with the TDHF results that are fitted by ujq ^ 65.5^"^/^ MeV. The evolution of the 
transition amplitude with A can be obtained from the energy weighted sum rule (EWSR) 
for quadrupole vibrations which reads [21] 



where {|a)} is an eigenbasis of H. In the last line of Eq. 30, we used the approximation of 
a constant density and a sharp surface which gives (r^) = j^R^ with R ki 1.2A^/^ fm. If 
all the strength is located at the GQR energy, which is a somewhat good approximation for 




Y.{E^-E,)\{a\QQm^ 



a 




14.3A^/=* MeV.fm^ 



(30) 
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FIG. 5: Evolution of (a) the GQR energy and (b) the transition ampHtude as a function of the 
number of nucleons in tin isotopes from the TDHF hnear response (circles). The line in panel 
(a) represents a A~^/^ fit of the TDHF results, whereas the line in panel (b) is obtained from 
considerations on the GQR energy weighted sum rule (see text). 

heavy nuclei [2] , then the EWSR reduces to 

S'Q^UQql^65.5A-^/YQ- (31) 
Equations (30) and (31) then lead to 

qq ^ 0.466A fm^. (32) 

This linear dependence is plotted in Fig. 5(b) and reproduces well the TDHF results. 

Let us now consider a dipole boost on these nuclei with a boost velocity ko — 0.01 fm~^. 
This value is small enough to generate a linear response of the dipole moment and a quadratic 
response of the induced quadrupole vibration in all considered isotopes. The GDR energy cud 
is shown as a function of the number of nucleons in Fig. 6(a). It is compatible with the 
dependence expected in heavy nuclei [2]. A fit of the TDHF results gives 

ojd ~ 76A-^/^ MeV. (33) 

Similarly to the quadrupole case, the dependence of the transition probability can be 



16 



100 110 120 130 140 

A 



FIG. 6: Evolution of (a) the GDR energy and (b) the transition probabihty as a function of the 
number of nucleons in tin isotopes from TDHF hnear response (circles). The line in panel (a) 
represents a A"^/^ fit of the TDHF results. The lines in panel (b) is obtained from considerations 
on the GDR energy weighted sum rule (see text) with an enhancement factor of the Thomas- 
Reiche-Kuhn sum rule n = 0.25 (dashed line) and n = 0.183 (solid line). 



obtained from the dipole EWSR 



^i-^d + 'Of (34) 



where k is the enhancement factor of the Thomas-Reiche-Kuhn (TKR) sum rule. Assuming 
all the strength in the GDR, i.e., qjj — S\,/ijJd, we get from Eqs. (33) and (34) 

2 ,^ , NZ 

In nuclear matter, the enhancement factor of the TKR sum rule \s n = 0.25 with the 
SLy4 parametrization [44]. This value clearly overestimates the transition probabilities (see 
dashed line in Fig. 6(b)), though the qualitative trend is in good agreement with the TDHF 
results. It is possible to compute k in finite nuclei using [50] 

'^ = ^^[^1(2 + xi) +h{2 + X2)] j drV„(r)pp(r) (36) 

where the U and Xi are the usual Skyrme parameters. In the considered tin isotopes, k is 
almost constant within the range 0.181—0.186 with no particular isospin or mass dependence. 
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This leads to a better agreement with the transition probabihties obtained with TDHF (see 
sohd hne in Fig. 6(b)), though a shght overestimation remains. The latter could be attributed 
to the fragmentation of the isovector dipole response. In this case, the dipole response reads 



In our calculations, the GDR properties {cud and qo) are extracted from the first minimum 
of the dipole response, which obeys to 



The last inequality implies that the TKR sum rule allows one only to compute the upper 
limit of qj) in our model. 

Finally, we investigate the coupling between the quadrupole and dipole vibrations from 
the quadratic response. We have shown in Sec. II A that, in the presence of a nonzero matrix 
element vq of the residual interaction coupling the state \D) to the state \DQ), a dipole 
boost is expected to generate an oscillation of the quadrupole moment. Before studying the 
evolution of vq along the tin isotopic chain, it is mandatory to get a deeper insight into the 
mechanism responsible for this induced quadrupole excitation. 

In a macroscopic approach, the isovector GDR is interpreted by a combination of 
the Steinwedel- Jensen model in which the total density is kept unchanged [51] and the 
Goldhaber-Teller model where proton and neutron fluids are incompressible [4]. It is obvi- 
ous that the Steinwedel- Jensen model does not affect the quadrupole moment because any 
modification of the density of one isospin specie is exactly compensated by the other in 
every point of space. In the Goldhaber-Teller model, however, a displacement of the proton 
and neutron spheres in the opposite direction is considered. It produces a dipole moment 
Qd = where X is the distance between their centers. This displacement also induces 

a prolate shape with a quadrupole moment quadratic in X . Indeed, assuming a displace- 
ment of a proton (resp. neutron) homogeneous sphere of density Zp^/A (resp. Npo/A) by 
Xp — —XN/A (resp. X^ — XZ/A) produces a quadrupole moment 



(37) 




(38) 



Qq oc ZXl + NXl 




(39) 



Using Eqs.(6) and (35), one gets Qq oc NZ/A^/^. 
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FIG. 7: Evolution of the coupling with mass number. The matrix clement of the residual interaction 
is plotted as a function of 1/A. The line shows a linear fit of the TDHF results. 

Finally, together with Eqs. (7), (32) and (35), the evolution of the coupling simply reads 
vq oc 1/A. This is, indeed, in agreement with the TDHF results shown in Fig. 7. It 
is interesting to note that, in this simple approach, the coupling does not depend on the 
isospin of the nuclei, but only on their total number of nucleons. In fact, the decrease of the 
absolute strength of the coupling with the number of nucleons is attributed to the fact that 
these couplings are mediated by the surface [16]. One then expects less anharmonicities in 
heavy nuclei. 

V. CONCLUSIONS 

We have shown that the residual interaction is responsible for anharmonicities in nuclear 
vibrations using three different analyses of time evolutions of multipole moments. We inves- 
tigated the coupling between one- and two-phonon states using a 3-dimensional TDHF code 
with a full Skyrme energy density functional. In particular, the excitation of a GDR couples 
to a GQR built on top of it, inducing a quadratic response of the quadrupole moment. The 
same coupling is responsible for the change of the GDR energy in static deformed states. 
The latter could be investigated using deformed RPA codes. As a consequence, the dipole 
frequency is modulated in case of dynamical deformation, e.g., induced by a quadrupole 
boost. This last property, associated with a Fourier analysis, might be used to investigate 
couplings when more than one mode is excited with the same quantum numbers. We finally 
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investigated these couplings with the quadratic response in several spherical tin isotopes. As 
a result, no dependence with isospin was found whereas an overall decrease of the coupling 
is obtained with increasing mass, showing that the couplings are mediated by the surface. 
These observations are interpreted within the Goldhabcr- Teller macroscopic model. These 
results indicate that no anharmonicity enhancement is expected in the standard giant reso- 
nances even for very exotic nuclei. However, couplings with exotic modes such as the pygmy 
dipole resonance should be investigated with the present method. The role of pairing and 
static deformation should be considered as well. 
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